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Abstract 

Real-space renormalization group maps, e.g., the majority rule transformation, map 
Ising-type models to Ising-type models on a coarser lattice. We show that each coefficient 
in the renormalized Hamiltonian in the lattice-gas variables depends on only a finite num- 
ber of values of the renormalized Hamiltonian. We introduce a method which computes 
the values of the renormalized Hamiltonian with high accuracy and so computes the co- 
efficients in the lattice-gas variables with high accuracy. For the critical nearest neighbor 
Ising model on the square lattice with the majority rule transformation, we compute over 
1,000 different coefficients in the lattice-gas variable representation of the renormalized 
Hamiltonian and study the decay of these coefficients. We find that they decay exponen- 
tially in some sense but with a slow decay rate. We also show that the coefficients in the 
spin variables are sensitive to the truncation method used to compute them. 



1 



1 Introduction 



Real-space renormalization group transformations were introduced to study critical behavior 
in Ising-type models. There has been extensive numerical study of these transformations, and 
there is a rich picture of how they are believed to behave. However, there are essentially no 
mathematical results on these transformations. The usual definition of these transformations 
is only formal since it involves an infinite- volume limit which must be proved to exist. The 
mathematical problem is to show that these renormalization group maps are rigorously defined 
in a neighborhood of the critical point, and to use them to study the system in a neighborhood 
of the critical point. This is a difficult problem and the amount of rigorous progress that has 
been made is embarrassing. Starting with the critical nearest neighbor Hamiltonian, the first 
step of the renormalization group transformation has been proved to be defined for a few specific 
lattices and transformations (TJ E]. The existence of the transformation well inside the high- 
temperature phase has been proved by rigorous expansion methods [2} El E] . It is possible to 
construct examples of transformations for which the renormalized Hamiltonian can be proved 
to be non-Gibbsian, including examples which start from the critical nearest neighbor Ising 

model manu. 

Even if we start with a finite-range Hamiltonian, after just one step of the renormalization 
group transformation the renormalized Hamiltonian will be infinite range and have infinitely 
many different terms. The conventional wisdom is that they should decay both as the number 
of sites involved grows and as the distance between these sites grows, so that the renormalized 
Hamiltonian may be well approximated by a finite number of terms. In some sense, this 
property is the raison d'etre of the renormalization group. It should allow one to study critical 
phenomena, which are inherently multiscale and so impossible to approximate well by a finite 
sets of terms, by studying a map of Hamiltonians which can be well approximated by a finite 
number of terms. 

Swendsen showed that one can compute the linearization of the renormalization group trans- 
formation about the fixed point from correlation functions that involve the original spins and 
the block spins [13]. His method allows one to avoid computing any renormalized Hamiltonians. 
From the point of view of using the renormalization group to calculate the critical exponents, 
this was a tremendous advance and was used in a large number of subsequent Monte Carlo 
studies of the renormalization group. From the point of view of trying to learn more about 
the renormalized Hamiltonians and the fixed point of the transformation, it had the unfortu- 
nate side effect that many of these Monte Carlo studies did not compute any renormalized 
Hamiltonians. In recent years there have been more studies that compute the renormalized 
Hamiltonian. In particular the the Brandt-Ron representation introduced in [JJ and studied 
further in [9j [101 El D2] is similar to the method we use in this paper. 

The goals of this paper are to give a highly accurate method for computing the renormalized 
Hamiltonian which works in the lattice-gas representation and to use it to test the conventional 
wisdom that the renormalization group transformation is well approximated by a finite number 
of terms. Our numerical calculations are done for the critical nearest neighbor Ising model on 
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the square lattice, and we only consider the renormalized Hamiltonian obtained by a single 
application of the majority rule transformation using 2 by 2 blocks. However, our approach 
is quite general and can be applied to other dimensions, lattices and choices of the real-space 
renormalization group transformation. 

One of the key tenets of the renormalization group is that if we fix a block-spin configuration 
and study the original system subject to the constraint imposed by the block spins, then this 
constrained system is in a high-temperature phase even if the unconstrained system is at its 
critical point. As an extreme case consider the block-spin configuration of all +l's with the 
majority rule transformation. The effect of this constraint on the original Ising system is similar 
to imposing a positive magnetic field, and the constrained system should have a relatively short 
correlation length. Our computational method for the renormalization group transformation 
takes advantage of this property. 

In the next section we review the definition of real-space renormalization group transforma- 
tions. In section three we explain our method for computing the renormalized Hamiltonian in 
the lattice-gas representation. Some of the details are postponed to section five. We use this 
method to study the decay of the terms in the renormalized Hamiltonian. In section four we 
consider how to compute the renormalized Hamiltonian in the more standard spin variables. 
There are multiple ways to do this, and we will see that the computed value of an individual 
coupling coefficient in the renormalized Hamiltonian varies considerable with the method used. 
We also study the decay of the renormalized Hamiltonian in the spin variables. Section five 
provides further detail for our method for computing the renormalized Hamiltonian. We con- 
sider the various sources of error in our computations in section six, and offer some conclusions 
in section seven. 

The significant dependence of coefficients in the renormalized Hamiltonian on the truncation 
method used has been seen before. In particular, Ron and Swendsen observed a change of several 
percent in the nearest neighbor coupling when the number of couplings kept was changed from 
six to twelve [9]. In [TO] they wrote "Even though the individual multispin interactions usually 
have smaller coupling constants than two-spin interactions, the fact that they are very numerous 
can lead to multispin interactions dominating the effects of two-spin interactions." Truncating 
the space of Hamiltonians implies that the linearization of the renormalization group map about 
the fixed point is also truncated. An early, interesting paper on the effect of this truncation is 

2 Real-space renormalization group transformations 

In this section we quickly review the definition of real-space renormalization group transforma- 
tions. We refer the reader to [16] for more detail. 

Consider an Ising- type model in which the spins take on only the values ±1. The lattice 
is divided into blocks and each block is assigned a new spin variable called a block spin. The 
example of the square lattice with 2 by 2 blocks is shown in figured! We consider transforma- 
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tions in which the block spins also take on only the values ±1. The transformation is specified 
by a kernel T(a,a). Here a denotes the original spins and a the block spins. The kernel is 
required to satisfy 



5> 



<7,<t) = 1 (1) 



for all original spin configurations a. The renormalized Hamiltonian H(a) is formally defined 
by 



Y,T(a,a)e- H ^ (2) 



(Note that the inverse temperature /3 has been absorbed into the Hamiltonians in the above 
equation.) This is only a formal definition since we must first restrict to a finite volume in 
order to make sense of this equation. Proving that the finite-volume definition of H has an 
infinite-volume limit is essentially an open problem. The condition (JTJ implies that 



so that the free energy of the original model can be recovered from the renormalized Hamilto- 
nian. This property allows one to study the critical behavior of the system by studying iterations 
of the renormalization group map. In particular, the critical exponents may be related to the 
eigenvalues of the linearization of the map about its fixed point. 

One widely studied family of kernels is the family of majority rule transformations. If there 
are an odd number of spins in every block, then T(a, a) = 1 if the majority of the spins in 
each block agree with the block spin and T(a, a) = otherwise. If there are an even number 
of spins in every block, then we let T(a, a) be the product over the blocks B of 

{1 if a B EieB °i > 
if ®b EieB °i < (3) 
1/2 if E^ = 

where <tb denotes the block spin for block B. 

The general approach presented in this paper applies to all these renormalization group 
maps. The numerical calculations that we will present are for the critical nearest neighbor 
Ising model on the square lattice with the majority rule renormalization group map with two 
by two blocks. 



3 Renormalized Hamiltonian in the lattice-gas variables 

Real-space renormalization group calculations are usually done using the spin variables <Tj = ±1. 
Our method is based on what are sometimes called the lattice-gas variables = (1 — <7j)/2 
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Figure 1: A renormalization group blocking for the square lattice. The original lattice 
sites are the open circles, the sites after the renormalization transformation are the 
solid circles. 



which take on the values 0, 1. Note that we have made the convention that a spin value of +1 
corresponds to a lattice gas value of 0. Throughout this paper we will use ex's for spin variables 
taking on the values ±1, and n's for lattice-gas variables taking on the values 0, 1. We indicate 
renormalized spins or variables with a bar over them, e.g., <7j, n^. We use a to denote the entire 
spin configuration {(Tj}. Likewise, n, a and n denote the corresponding collections of variables. 

In this section we work entirely in the lattice-gas variables, both for the original Hamiltonian 
and the renormalized Hamiltonian. We write the renormalized Hamiltonian as 

H(n) = Y,c(Y)n(Y) (4) 

Y 

where the sum is over all finite subsets including the empty set, and 

n{Y) = \[n l 



5 



Consider the block- variable configuration of all O's. Our method for computing the renor- 
malized Hamiltonian uses only block- variable configurations which differ from this configuration 
at a finite number of sites. For a finite subset X, let n x denote the block- variable configuration 
with all block variables in X equal to 1 and the rest equal to 0. Then eq. (J2]) says 

exp(-H(n x )) = J2 T(n x ,n) e~ H{n) 

n 

Note that n (X) = except when X = 0. So H(n <l> ) = c(0). In particular, c(0) will grow 
as the size of the finite volume. The other coefficients c(Y) should have finite limits in the 
infinite- volume limit. We define f(X) by 

/(X) = H(n x ) - H(n®) 

Then f(X) should have a finite limit in the infinite- volume limit, and it should be related to 
the infinite-volume c(X) by 

nx)= j2 < Y ) ( 5 ) 

Y-S^YCX 

The system of equations (EJ can be explicitly solved for the c(Y). We claim that the solution 
for X ^ is 

cW= £ (-l) |XHy| /(F) (6) 
This is a standard inversion trick. To verify ([6]), define c(X) by (E]). Then for a given X ^ 0, 

E < Y ) = E E (-i) |y|HZ| /(^) 

= E E (- 1 ) |y|HZ| ( 7 ) 

Z-M^ZCX Y-.ZCYCX 

The sum over Y is 1 if X = Z. If Z is a proper subset of X, we claim this sum is 0. To see 
this: 

e (-i) |yHzi = e (-!) iwr| = n (i+(-i))=o 

Thus (J71) collapses to f(X). 

The important feature of eq. ([6]) is that the coefficient c(X) only depends on a finite number 
of free energies f(Y), specifically those with Y C X. As we will see, these free energies can be 
computed extremely accurately. So individual coefficients c(X) in the lattice-gas variables can 
be computed extremely accurately. Moreover, this computation does not depend on how many 
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terms we decide to keep in the renormalized Hamiltonian. If we increase the number of terms 
we keep, then the coefficients we have already computed will not change. 

We have carried out numerical calculations of a large number of the coefficients in the 
lattice-gas representation for the critical nearest neighbor Ising model on the square lattice 
with the majority rule renormalization group map with two by two blocks. We need a criterion 
for deciding for which Y we will compute c(Y). We assume the coefficients will decay as the 
number of sites in Y grows and as the distance between these sites grows. So we need a measure 
of the size of a set Y. There is no canonical way to define this size. We use the following ad 
hoc quantity. If Y = {yi, 1/2, ■ ■ ■ , y n }, then we define 

n 

sor) = Y,\\vi-c\\l (8) 
1=1 

where c is the center of mass: 

1 - 

c = - V Vi 
n / — ' 

and || 1 1 2 is the usual Euclidean distance in the plane. Note that we do not take a square root 
in flHD. 

We claim that if X' C X then S(X') < S(X). To prove this it suffices to prove the case 
that X' has one less site than X. Let X be X\,X2, ■ ■ ■ ,x n and X' be X\, X2, ■ • • , x n -i- Let c be 
the center of mass of X and d the center of mass of X'. So 

n— 1 

S(X') = J2\\^-c'\\l 

»=i 

The center of mass has the property that it minimizes the above sum. So 

n— 1 n 

s(x') < in - c\\i < \\ x i - °wi = s ( x ) 

i=l i=l 

We fix a cutoff C > 0, and compute c(Y) for all Y with S(Y) < C . We only need to 
compute it for one Y from each translation class, and so there are a finite number of such 
Y's. The bulk of the computation is computing the free energies f{X) for X with S(X) < C. 
Using eq. (J6]) to find the c(Y) requires comparatively little computation. The property that 
I'cl^ S(X') < S(X), implies that the collection of X for which we must compute f{X) 
is just the collection of X with S(X) < C. 

To study how fast the coefficients c(Y) decay, we take one coefficient from each translation 
class that we have computed and order them so that \c(Y)\ is decreasing, i.e., |c(F„)| > |c(F„, + i)|. 
We then plot |c(F n )| as a function of n. This is the bottom curve in figure El Note that the 



7 



16 " Ub 2000 4000 6000 8000 10000 

n 

Figure 2: The coefficients c(Y n ) are ordered so |c(l^)| decreases. The bottom curve 
is |c(Y^)| vs. n, and the top curve is the tail J2iL n l c (^«)l vs - n - 

vertical axis uses a logarithmic scale. The second quantity plotted (the top curve in the figure) 
is 

N 
i=n 

as a function of n, where N is the total number of Y for which we compute the coefficients. 
The two lines shown are given by c2~ n / 850 for two different values of c. The two curves in the 
figure depend on the function S(Y) we use to measure the size of sets and the cutoff we use for 
this function. However, whatever function and cutoff we use, the resulting curve will be a lower 
bound on the curve that would result from computing all the coefficients c(Y). In particular, 
we can make the following observations. The lower curve crosses the horizontal lines at 10~ 2 , 
10 -3 and 10 -4 at 131,1223 and 4023, respectively. Hence there are at least 131 translation 
classes with a coefficient bigger than 10~ 2 , at least 1223 with a coefficient bigger than 10~ 3 , 
and at least 4023 with a coefficient bigger than 10~ 4 . 

In the preceding discussion we used one coefficient from each translation class. In addition 
to the translation symmetry the model is also symmetric under rotations by 90 degrees and 
relections in lattice axes. More precisely, the additional symmetry is the dihedral group of order 
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8. We have repeated the previous study of the decay of the coefficients taking into account the 
dihedral group symmetry as well as the translational symmetry by taking only one term from 
the above list from each dihedral group symmetry class. The main effect is that the scale on 
the horizontal axis is reduced by a factor of 8. This is not surprising since for most subsets, 
rotations and reflections will generate eight different subsets. 

From a mathematical perspective, one would like to show that the renormalized Hamiltonian 
exists in some Banach space. One choice of norm for the Banach space would be 

E i c ( y )i 

Y:0eY 

One would like to approximate the Hamiltonian by a Hamiltonian with a finite number of 
terms. So it is important to see how fast the above sum converges as we include more terms 
in the Hamiltonian. This is similar to the second plot in figure EJ Note that in this norm each 
translation class appears \Y\ times. So the second plot in figure [2] is in some sense a lower 
bound on the decay for the Hamiltonian. Another choice of the norm would be 

E HY)K {Y) 

Y-.oeY 

for some measure n(Y) > of the size of Y . For this norm the convergence would be even 
slower than that seen in the figure. 

It is worth noting that norms defined using the lattice-gas representation of the Hamiltonian 
are in general stronger than norms that use the spin variable representation [3]. For example, 
if we can write the Hamiltonian in the lattice-gas representation (J4]) with 

E icooi < oo 

Y30 

then it is straightforward to show that the Hamiltonian can be represented in the spin variable 
representation (Q with 

Eicni<°o 

Y30 

4 Renormalized Hamiltonian in the spin variables 

In the previous section we saw that in the lattice-gas variables there is a natural way to compute 
the coefficients c(Y) in the expansion (j3J) for H. In this section we consider the renormalized 
Hamiltonian in the spin variables: 

H(a) = E d(Y)a(Y) (9) 

Y 
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with a(Y) = YlieY The sum over Y is over all finite subsets. 

We can use = (1 — ^)/2 to express the spin coefficients d{Y) in (J9]) in terms of the 
lattice-gas coefficients c(Y) in (@|. 



= ^c(IHI)^c(I)2-WH(1-^ 

= E c W 2 ~' x| E (-^'^cn 



X Y:YCX 

= E-( y )(- 1 ) |y| E ^)2" w 

y x : ycx 

= E d ( y )^( y ) 

y 

where the spin coefficients d(Y) are given by 

<n = (-i) |y| E ^)2- |x| (io) 

x : ycx 

The problem is that to compute the spin coefficient d(Y) we need the lattice-gas coefficients 
c(X) for infinitely many X, and so we need the free energies f(X) for infinitely many X's. So 
we must introduce some sort of approximation. 

Let Yqo be a collection of finite subsets of the renormalized lattice such that one set from 
each translation class is contained in Y™. We can rewrite (19]) as 



H{a)= E d(Y)Y,HY + t) 



ye¥ 



where the sum over t is over the translations for the renormalized lattice. Here Y + t denotes 
{i + t : % G Y}. 

Now let Y be a finite subcollection of We want to compute an approximation to the 
above of the form 



iT(a)^E rf ( F )E (T ( r + t ) 



YeY t 

We will consider two methods which we will refer to as the "partially exact" method and the 
"uniformly close" method. 

For the partially exact method, we begin by noting that we can write H(n) as 

H(n)= E c(Y)J2n(Y + t) 

yeYoo t 

The approximation is simply to truncate this sum by restricting Y to those in Y: 

yeY t 
10 



The c(Y) are exact. As we saw in the last section we can compute them from ([6]) by computing 
the free energies f{X) for XgY. We then convert this Hamiltonian to the spin variables with 
no approximation. The result is that the approximation to H(a) is 



(11) 



ye¥ t 



with 



d(Y) = (-1) 



Y 



£ c(X) 2-1 



x 



(12) 



X:YdX,X+t&{ 



where the notation X + 1 G Y means some translation of X (possibly X itself) is in Y. Thus 
this method is equivalent to truncating the exact formula (TTUT) by restricting the sum over X 
to sets in Y and their translates. In the lattice-gas variables our approximation to H agrees 
with the true H for all n x such that X G Y. The change from lattice gas to spin variables did 
not involve any approximation, so our approximation to H in the spin variables agrees exactly 
with the true H for all configurations a Y for Y G Y. This is the reason for calling this method 
"partially exact." It is exact for some of the block-spin configurations. 

For the uniformly close method let X be another finite collection of finite subsets which 
contains at most one set from each translation class. We compute the free energies f{X) for 
X G X, i.e., we compute H{a x ). We define the error of a set of coefficients {d(Y) : Y G Y} to 



where a is the spin configuration which is —1 on X and +1 on all other sites. We then choose 
the coefficients d(Y) to minimize the above error. This is a standard linear programming prob- 
lem which we solve by the simplex algorithm. We call this the uniformly close approximation 
since we have a uniform bound on the difference between our approximation and the exact H 
for the block-spin configurations a x for X G X. (For other X we cannot say anything about 
how well the approximation does.) If X = Y, then the partially exact approximation makes the 
above error zero. We only use the uniformly close approximation for X which are larger than 



We take the following point of view. We think of Y as being fixed. It determines a finite- 
dimensional space of Hamiltonians that we use to approximate the renormalized Hamiltonian. 
We then think of the collection X as being variable. A larger X means we "know" more free 
energies and so have more information to use in computing the approximation. In our studies 
we will take the collection Y to be all the subsets Y with s(Y) < Cfj for some cutoff Cg, and 
X to be all the X with s(X) < Cf for some cutoff Cf > Cg 

When we worked in the lattice-gas variables the computation of the coefficients c(X) was 
unambiguous. The computation of the values of H(n) requires some approximations, but as 



be 




Y. 
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we will see in section [5] these approximation are well behaved and introduce small errors. The 
computation of the c(X) from the H(n) does not require any approximation or truncation. In 
the spin variable representation we now have multiple ways to compute the coefficients d(X) 
depending on whether we use the partially exact or uniformly close methods and on the choices 
of X and Y. We restrict our study of the spin variable coefficients to studying how these 
choices affect the values of individual coefficients. We focus our attention on three particular 
coefficients: the nearest neighbor, the next nearest neighbor and the plaquette. These refer 
to the coefficients of a^aj with \i — j\ = 1, of a^aj with \i — j\ = v2, and of a^aja^ai where 

k, I are the corners of a unit square. As in the previous section, our numerical calculations 
are for the critical nearest neighbor Ising model on the square lattice with the majority rule 
renormalization group map with two by two blocks. 

For the partially exact method we have one parameter - the cutoffs Cg and Cf are equal 
and correspond to the cutoff C of section [31 So we can plot individual coefficients as a function 
of Cf. For the uniformly close method we have two parameters: the cutoff C# determines the 
finite-dimensional subspace used to approximate the renormalized Hamiltonian and the cutoff 
Cf determines the number of H{n) values we use. We plot the coefficients in this case as a 
function of Cf for several different choices of Cg. The results are shown in figures |3"]H^ and[5j 
The variations seen in the three coefficients are roughly comparable in size. Note that while 
the ranges of the vertical axes vary in the three figure, the scales for the vertical axes are all 
the same. The variations in the coefficients are on the order of several thousandths. So even 
for these relatively large coefficients, it is difficult to determine the value of the coefficient to 
better than a few percent. For smaller coefficients the variations are somewhat smaller, but as 
a fraction of the coefficient they are typically larger. 

These three coefficients (along with many others) have of course been computed before. 
Two early references are 0, EE]. The point of our study is not the values of these coefficients 
but rather the variation in their values as one varies the method used to compute them. 
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Figure 3: The dependence of the nearest neighbor coefficient in the spin variable 
representation of the renormalized Hamiltonian on the computation method. The 
solid curve is the partially exact method. The dashed curves are the uniformly close 
method with four different choices of Cg 
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Figure 4: The next nearest neighbor coefficient in the spin variable representation of 
the renormalized Hamiltonian 
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Figure 5: The plaquette coefficient in the spin variable 
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5 Computing the free energy 



Fix a block-spin configuration a. We want to compute the free energy H{a) of the constrained 
partition function 

exp(-jf(*)) = X;T(*,<7)e-*« 

(j 

Initially we work with the spin variables, but later in this section we will switch to the lattice- 
gas variables. We only need to do this computation for configurations a which are +1 except 
on a finite set. Even when the original system is at the critical point, these constrained systems 
have relatively short correlation lengths. This is where the real power of the renormalization 
group becomes apparent. In particular, finite- volume effects in the above computation decay 
very quickly as the volume increases. 

Before we explain our method for the computation, we first indicate how H(a) can be 
computed by a Monte Carlo calculation. (We have done such a simulation as a check on the 
method we describe later.) Fix a relatively small set of block spins, X. Let V be a finite 
volume of block spins containing X which is large enough that the boundary of V is far from 
X. We include only the factors in the renormalization group kernel corresponding to the blocks 
in V \ B. For these blocks we take the block spins to be +1. We then run a Monte Carlo 
simulation of the Ising system with this kernel outside of X. When we sample the simulation 
we compute the block-spin configuration on X. This allows us to compute the relative weights 
of the possible block-spin configurations on X. From these weights we can then compute the 
H{a) for a which are —1 only on a subset of X. 

We now turn to our method for computing H(a). It does not involve Monte Carlo meth- 
ods, and it is much more accurate than the Monte Carlo approach described in the previous 
paragraph. Everything in the following depends on a, but we will not make this dependence 
explicit. Imagine that we have summed over the spins one block at a time in such a way that 
we have reached the state in figure [6j Open circles indicate sites in the original lattice for which 
we have already summed over the spin, and blue (gray) circles represent sites for which we have 
not. (Red (solid) circles indicate the block spins which are fixed throughout this computation.) 
The result of this partial computation of the free energy is a function of the spins in the original 
lattice with shaded circles. In fact, it only depends on those that are nearest neighbors of a 
spin with an open circle. We will refer to these spins as boundary spins. The quantity we have 
computed so far is positive, and we write it in the form 

exp(^a(X)cx(X) + AH) Y[t B ,(a,a) 

X B> 

where X is summed over finite subsets of the boundary spins. AH denotes the terms in the 
Hamiltonian that only involve spins with shaded circles. (These terms have not yet entered the 
computation.) The product over B' is over the blocks containing shaded circles, and tB'(&,cr) 
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Figure 6: We compute H{a) by summing out the original spins one block at a time. 
The open circles are spins that have been summed over, while the blue (gray) circles 
are spins that have yet to be summed over. The red (black) circles are the fixed block 
spins. 



is the factor in the renormalization group kernel for block B'. The next step is to sum over the 
four spins in the block B and take the logarithm of the result: 



In 



^expQTa(Xy(X) + Atf)n 



X 



13' 



The sum over a B denotes a sum over the spins <7j with i G B. Terms a(X)a(X) for which 
X n B = pass through this computation trivially. So do the terms in AH which do not 
involve a spin in the block B and the factors t B ' for B' ^ B. So the computation that we must 
actually do is 



In 



^exp( a{X)a{X) + h)t B (a,a) 

o- B X:XnB^9) 
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where h contains the terms in H that only depend on spins with shaded circles and depend on 
at least one spin in B. 

To do this computation numerically, we must introduce a truncation. We fix a finite subset 
D of the boundary sites centered near B. We then restrict the sum over X to X C D. We 
need to write the result of the truncated computation in the form 



In 



^exp( V~ a(X)a(X) + h) t B (a,a) 

o-b X:XnB^9,XcD 



Y 



The left side only depends on spins in D' = D \ B, so the sum on the right may be restricted 
to Y C D' . If we define F(a) to be the left side of this equation, then the coefficients are given 
by 



a'(Y) 



,-|D 



cr 



The amount of computation required grows quite rapidly as D grows for three reasons. First, 
the number of X with X C D grows as 2' D l Second, the sum over gd 1 in the above also grows 
as 2' Z5 L Third, the number of Y also grows as 2' jD L We have found that a(X) decays quickly 
as the number of sites in X grows. So we can make a further truncation by only keeping terms 
a(X) with \X\ less than some specified cutoff. (\X\ denotes the number of sites in X.) This 
greatly reduces the growth of the computation with D from the first and third effects. But we 
are still left with the second effect. 

We can eliminate the second effect by working in the lattice-gas variables. We replace 
Zx a(X)a(X) by £ x b(X)n(X). Define 



F(n) = In 



^exp( b{X)n{X) + h)t B \ 

n B X:XnB^9,XcD 



n,n) 



(13) 



We need to compute the coefficients in 

F(n) = J2b'(X)n(X) 



x 



As we saw in section [31 they are given by 



b'(X)= £ (-l)W-^F(^) 

Y-S^YCX 



(14) 



where n Y is the configuration that is 1 on Y and off of it. So to compute b'(X) we only need 
to compute F{n Y ) for Y C X. 
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In this approach using the lattice-gas variables we can forget about the set D entirely. 
Instead we specify a finite collection B of subsets of the boundary spins with the property that 
they intersect B. We then make the approximation 



(15) 



X:XnB^<A 



We use (tHJ) to compute b'(X). It will be nonzero only for X C D' . Before we sum over the 
next block of spins, we need to truncate b' (X)n(X) . We keep only the terms such that X 
is in B + 1 where t is the translation that takes the block we just summed over to the block we 
are summing over next, and B + t denotes the collection of sets of the form X + t for X G B. 

We take the finite collection B to be all X which intersect B and satisfy S(X) < Cb where 
S(X) is some size function and Cb is some cutoff. We use the size function given by flS} that 
we used for choosing the block-spin sets. In our calculations we take Cb = 260 which leads to 
10, 763 sets in the collection B. We discuss the effect of Cb on the error in the following section. 

The above discussion took place in an infinite volume. The region shown in figure [6] is a 
finite piece of the infinite volume. In practice we can only sum over the spins in a finite number 
of blocks. The block-spin configuration n is of the form n Y for a finite set Y. We carry out the 
computation in a finite volume which is chosen so that the distance from Y to the boundary 
of the finite volume is sufficiently large. We will study how large the finite volume should be 
in the next section. 

6 Errors 

In this section we study the sources of error in our computations of the f{X). There are two: 
the use of a finite volume to compute infinite-volume quantities and the truncation determined 
by the cutoff Cb in section [51 

We choose the finite volume in which we carry out our calculation as follows. The block-spin 
configurations n that we consider are of the form n Y for finite sets Y . We take these sets Y 
to be centered near the origin and take the finite volume to be a square centered at the origin. 
The square contains the blocks with centers at (2i, 2j) with — L < i < L, —L < j < L. So the 
infinite- volume limit is obtained by taking L — > oo. 

To study the finite-volume error in our calculation we do the following. The free energy 
f(Y) depends on L, so we denote it by fiiY). As a measure of the finite-volume error we use 



where the sum is over one element of each translation class with s(Y) < 210, and N is the 
number of terms in the sum. In this study of the finite- volume error we take Cb = 30. This 
is much smaller than the cutoff we used for the main calculations, but the behavior of the 




(16) 



Y 
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Figure 7: L is a measure of the size of the finite volume. The quantity plotted is the 
average change in f(X) when L is decreased by 1. (See eq. (fl6l).) 



finite-volume error with L is the same whether we look at this smaller set of coefficients or the 
larger set. 

This average difference as a function of L is shown in figure [7J The vertical scale is log- 
arithmic, so the approximately linear dependence seen for the smaller values of L indicates 
exponential decay of this difference with L. The line shown in the figure is of the form ce _L// °' 2 . 
Keeping in mind that L corresponds to numbers of blocks and the blocks are 2 by 2, the decay 
length of 0.2 corresponds to a decay length of 0.4 in units of lattice spacings. This very short 
decay length is a result of the block spin being +1 at all but a finite number of block sites. 
Beginning with L around 10 the difference is dominated by numerical errors. In our simulations 
we are very conservative and take L = 15. With this choice the finite-volume error is at the 
same level as the numerical error. 

The translational symmetry of the original model implies that f(X) is unchanged if we 
translate X. So we only need to compute f(X) for one element of each translation class. The 
model is also invariant under the dihedral group symmetry generated by rotations by mutiples 
of 7r/2 and reflections in the coordinate axes. However, our method of computing f(X) breaks 
the dihedral symmetry of the lattice, so our values of f{X) for X's from the same dihedral class 
are not exactly the same. The dihedral symmetry is only restored when we let L — > oo and 
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Figure 8: We plot the average of f(X) — f(X), as defined by eq. (ITT]) , and of 
c(X) — c(X), as defined by eq. (fTBl . as a function of the cutoff Cb- 



Cb — > oo. We have already seen that we can take L sufficiently large that the finite- volume 
error is reduced to the order of the numerical error. So we can use the breaking of the dihedral 
symmetry to study how the error depends on the cutoff Cg. 

For various choices of we compute f{X) for the same collection of X as in our main 
calculation. Let f{X) be the average of f(Y) over one Y from each translation class which is 
related to X by the dihedral symmetry. (The number of terms involved in this average ranges 
from 1 to 8.) The differences f{X) — f(X) are a measure of the amount of breaking of the 
dihedral symmetry and hence of the error in the computation from the cutoff Cb- We use the 
average 

^Ei/cn-zcni ( 17 ) 

Y 

to quantify the error. As before the sum is over one element of each translation class with 
s(Y) < 210, and N is the number of terms in the sum. This quantity is plotted in figure El as 
a function of Cb for the free energies f(Y). It is the higher set of points. For the coefficients 
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in the lattice-gas variables we define c(X) analogously, and study the average 



This quantity is the lower set of points in figure [HJ 
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Figure 9: The convergence of four different quantities as Cb — > oo. From top to 
bottom the four quantities are given by equations (fT9j) to (1221) . 



We also study the convergence as Cb — > oo in another way. Let f°°(Y) denote f(Y) for the 
largest value of Cb which we use, i.e., 260. We then consider 



(19) 



This is plotted as a function of Cb in figure [9] for the free energy f(Y). We also plot 

jj'EifM-rool (20) 
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As the figure shows, averaging over the dihedral group like this reduces the error somewhat. 
The figure also includes the analogous plots for the coefficients in the lattice-gas representation, 
i.e., of the quantities 

Lj2\c(Y)-c°°(Y)\ (21) 

Y 

and 

l£|g(y)-e~(y)| (22) 

7 Conclusions: 

We have shown that if we use lattice-gas variables, then the computation of the coefficients in 
the renormalized Hamiltonian only depends on a finite number of values of the renormalized 
Hamiltonian. So this computation does not depend on how we approximate the inherently 
infinite-dimensional renormalized Hamiltonian by a finite-dimensional approximation. We have 
also given a highly accurate method for computing the values of the renormalized Hamiltonian 
which takes advantage of the finite correlation length that results from the introduction of the 
renormalization group transformation. 

The renormalized Hamiltonian has infinitely many different terms but the conventional 
wisdom is that it may be well approximated by a finite number of terms. In particular, the 
magnitude of the coefficients should decay as the "size" of the set of lattice sites increases. 
We studied this for the nearest neighbor critical Ising model on the square lattice under one 
step of the majority rule renormalization group transformation. We computed a large number 
of coefficients in the lattice-gas variables, ordered them by decreasing magnitude and plotted 
them. We found that over several orders of magnitude the coefficients decayed exponentially 
with the number of terms, but the decay rate was slow. It takes about 850 additional terms to 
see the magnitude reduced by just a factor of 1/2. 

If we use the usual spin variables, there is no natural way to compute the coefficients of 
the renormalized Hamiltonian. We considered two methods of truncation. If we look at an 
individual coefficient, we see significant dependence on the method used and on the value of 
the cutoffs used to specify the truncations in these methods. Even with our computation 
of approximately 10, 000 values of the renormalized Hamiltonian, the uncertainty in the spin 
variable coefficients due to the different truncation methods is on the order of a percent for the 
largest coefficients and even larger as a percentage for some of the smaller coefficients. 

One might hope to prove theorems about these real-space renormalization group transfor- 
mations by defining a suitable Banach space of Hamiltonians and then doing a computer aided 
proof to show the transformation is defined in some open subset of the Banach space and there 
is a fixed point in this subset. Proving there is a fixed point would require constructing an 
approximation to the fixed point with a finite number of terms. Our numerical results suggest 
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that at best such an approach will require a huge number of terms in the finite approximation 
and at worst the number of terms needed will doom the approach to failure. 

Past numerical studies of the two dimensional Ising model using the renormalization group 
have produced fairly accurate values of the critical exponents using a relatively modest number 
of terms in the renormalized Hamiltonian. These studies use the spin variables, so their accuracy 
is surprising given the difficulty we have found in computing the coefficients in the renormalized 
Hamiltonian accurately. An interesting question is to understand this. 



References 

[1] A. Brandt and D. Ron, Renormalization multigrid (RMG): statistically optimal renor- 
malization group flow and coarse-to-fine Monte Carlo acceleration, J. Stat. Phys. 102, 
231-257 (2001). 

[2] R. B. Griffiths and P. A. Pearce, Mathematical properties of position-space 
renormalization-group transformations, J. Stat. Phys. 20, 499-545 (1979). 

[3] R. Gupta and R.Cordery, Monte Carlo renormalized Hamiltonian Phys. Lett. 105 A, 415- 
417 (1984). 

[4] R. B. Israel, Convexity in the theory of lattice gases (Princeton University Press, Princeton, 
1979). 

[5] R. B. Israel, Banach algebras and Kadanoff transformations, in: Random Fields (Es- 
ztergom, 1979), vol. II, J. Fritz, J. L. Lebowitz, D. Szasz, (eds.), (North-Holland, 1981). 



[6] I. A. Kashapov, Justification of the renormalization - group method, Theor. Math. Phys. 
42, 184-186 (1980). 

[7] T. Kennedy, Some rigorous results on majority rule renormalization group transformations 
near the critical point, J. Stat. Phys. 72, 15-37 (1993). 

[8] T. Kennedy and K. Haller, Absence of renormalization group pathologies near the critical 
temperature - two examples, J. Stat. Phys. 85, 607-637 (1996). 

[9] D. Ron and R. H. Swendsen, Calculation of effective Hamiltonians for renormalized or 
non-Hamiltonian systems, Phys. Rev E 63, 066128 (2001). 

[10] D. Ron and R. H. Swendsen, Importance of multispin couplings in renormalized Hamilto- 
nians Phys. Rev E 66, 056106 (2002). 



24 



D. Ron,R. H. Swendsen and A. Brandt, Inverse Monte Carlo renormalization group trans- 
formations for critical phenomena, Phys. Rev. Lett. 89, 275701 (2002). 

D. Ron,R. H. Swendsen and A. Brandt, Computer simulations at the fixed point using an 
inverse renormalization group transformation, Physica A 346, 387-399 (2005). 

R. Shankar, R. Gupta, and G. Murthy, Dealing with truncation in Monte Carlo 
renormalizat ion-group calculations, Phys. Rev. Lett. 55, 1812-1815 (1985). 

R. Swendsen, Monte Carlo renormalization group, Phys. Rev. Lett. 42, 859-861 (1979). 

R. Swendsen, Monte Carlo calculation of renormalized coupling parameters. I. d = 2 Ising 
model, Phys. Rev. B 30, 3866-3874 (1984). 

A. C. D. van Enter, R. Fernandez and A. D. Sokal, Regularity properties and pathologies of 
position-space renormalization - group transformations: scope and limitations of Gibbsian 
theory, J. Stat. Phys. 72, 879-1167 (1993). 

A. C. D. van Enter, Ill-Defined Block-Spin Transformations at Arbitrarily High Temper- 
atures, J. Stat. Phys. 83, 761-765 (1996). 



25 



